Dissipation control in cavity QED with oscillating mode-structures 
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We demonstrate how a time-dependent dissipative environment may be used as a tool for con- 
trolling the quantum state of a two-level atom. In our model system the frequency and coupling 
strength associated with microscopic reservoir modes are modulated, while the principal features of 
the reservoir structure remain fixed in time. Physically, this may be achieved by containing a static 
atom-cavity system inside an oscillating external bath. We show that it is possible to dynamically 
decouple the atom from its environment, despite the fact that the two remain resonant at all times. 
This can lead to Markovian dynamics, even for a strong atom-bath coupling, as the atomic decay be- 
comes inhibited into all but a few channels; the reservoir occupation spectrum consequently acquires 
a sideband structure, with peaks separated by the frequency of the environmental modulation. The 
reduction in the rate of spontaneous emission using this approach can be significantly greater than 
could be achieved with an oscillatory atom-bath detuning using the same parameters. 

PACS numbers: 42.50.Ct, 03.65.Yz 
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I. INTRODUCTION 

It is well known that the ultimate origin of dissipa- 
tion is the unavoidable interaction between a system and 
its surroundings. For quantum systems, a distinction is 
usually made between decay, during which energy leaks 
into the environment, and decoherence, whereby correla- 
tions established with the surroundings rapidly wash-out 
any coherence between certain elements of the system 
[l[ . Decoherence generally takes place on a much shorter 
timcscalc than decay, and this timescale is governed by 
the separation in phase-space of different components of 
the system density operator. However, for a single two- 
level atom, the timescale governing both decay and de- 
coherence is related to the rate of spontaneous emission, 
which simplifies matters and thus makes this an ideal 
quantum system in which to study dissipative effects. 

The idea of using the surroundings to control atomic 
decay has been pursued ever since the first studies using 
nearby conducting surfaces and cavities 0, d, 0, H[ and 
experimental control has now reached the point where an 
atom can be strongly-coupled to a cavity in which it is 
placed. This technological development has allowed for 
many useful applications of the atom-cavity interaction, 
including the generation of entanglement I6| , creation of 
number states of the electromagnetic field [3, HI and pho- 
tons on demand observation of quantum jumps of 
the electromagnetic field [ll[ and experimental tests of 
non-locality fl2| among many others [l3|. Nevertheless, 
interaction with the environment is traditionally seen as a 
negative feature since this is the root-cause of dissipation 
and there is currently a huge theoretical and experimen- 
tal effort aimed at controlling dissipative effects in cavity 
quantum optics. 

One way to achieve this goal is to limit the extent of the 
system-environment coupling, although ultimately this 
will always still prevail at some level. Another promising 
approach is to exploit certain symmetries in the system- 



environment coupling, so as to confine the dynamics to a 
(decoherence-free) subspace of the overall Hilbert space 
in which dissipative effects cancel out [HI, [H| • A third 
method, and the approach we shall pursue here, is to 
dynamically modify the system's surroundings so as to 
control the effects of the dissipative couplings rather than 
reducing their magnitude [l6|, EE d, O HH ■ The use of 
dynamically engineered reservoirs to control dissipation 
is appealing, since this technique is non-invasive and fur- 
thermore does not require individual components of the 
system to be addressed separately. 

To the lowest order of approximation, an atom's en- 
vironment consists simply of the electromagnetic field 
modes into which it can emit. We should therefore expect 
that changing some property of all of these modes should 
affect the process of atomic decay, and these effects are 
well documented (l8l . [2ll . |22| . It is less clear what will 
happen in the more subtle case where the properties of 
individual modes are altered, but in a controlled way, 
such that, for example, the combined effect is to keep 
all macroscopic reservoir structure static in time. It is 
tempting to think that this type of manipulation might 
only produce observable effects in strongly- coupled atom- 
cavity systems, since it is only the memory kernel for the 
atom-bath interaction which changes and is dynamically 
modified under this model. However, it has recently been 
shown (in Ref. [23j for a linear increase of all reservoir 
frequencies) that this form of reservoir manipulation can 
be used to control the rate of spontaneous emission for 
both strong and also weak atom-reservoir coupling. 

In the current paper we investigate the case of an oscil- 
latory manipulation of all reservoir mode frequencies and 
coupling strengths. This can be achieved by modifying 
the length of a cavity, and hence the mode frequencies, 
in an oscillatory way, (further details are given in sec- 
tion [V] of this article). From a practical perspective, an 
oscillating mode-structure has several advantages over a 
linear chirp of the reservoir frequencies: (i) piezoelectric 
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actuators used to modify the cavity length typically per- 
form favourably in an oscillatory regime; (ii) the required 
changes in the cavity length are periodic and do not grow 
large at long times, unlike in the linear-chirp case. 

The remainder of the article is arranged as follows. 
In section UH the mathematical model is introduced and 
some general observations are made regarding the appli- 
cability of dynamic environments to the control of deco- 
herence and decay. In section UTTl the detailed form of the 
microscopic reservoir structure is given for an oscillatory 
modulation of the reservoir mode frequencies and closed- 
form expressions for the atomic decay rate and emission 
spectrum are derived. These results are tested against 
numerical simulations in section IIVI for the specific case 
of a sinusoidal modulation of the bath mode frequencies. 
A particular physical realisation of the model and the 
potential size of observable effects are examined in detail 
in section [Vj It is found that our scheme can give rise to 
a surprisingly large reduction in the rate of dissipation 
and that the method studied here compares favourably 
with other, more basic reservoir manipulations studied 
previously (e.g. 0, [H, EI El ) ■ Finally, in section Efl 
we conclude our findings. 



II. GENERAL DYNAMIC ENVIRONMENTS 
FOR A TWO-LEVEL ATOM 

We consider a two-level atom with transition frequency 
luq and lower and upper states |0) and |1), coupled to 
a zero-temperature bath of electromagnetic field modes, 
which together constitute a reservoir for the atomic de- 
cay. The reservoir is engineered so that the individual 
bath mode frequencies uik(t) and also the coupling g k (t) 
between the atomic transition |0) «-> |1) and the k th mode 
of the radiation field are time- dependent. We initially 
consider a discrete bath. Without loss of generality, all 
couplings are chosen to be real, since any time-dependent 
phase in the coupling g k (t) can be transferred onto the 
time-dependence of the corresponding mode frequency 
tu k (t) and as usual, constant phase terms can be absorbed 
into the basis states. The Hamiltonian for the composite 
atom-reservoir system in the rotating-wave approxima- 
tion (with h = 1) is: 



H(t) = LU Q a+a- + Wfc(«) (b{b k + 1/2J 

k 



with raising and lowering operators er + = |1)(0 and 
(7~ = |0)(1|. Due to the oscillation of the cavity, indi- 
vidual mode-frequencies may pass through resonance on 
a timescale which is very rapid from the point-of-view of 
the atom and may even be much shorter than the bath 
correlation time. However, we note that in any practical 
setting, the rate of change of reservoir mode frequencies 
is never fast enough to create photons in the reservoir by 



itself (in the sense of refs. [2y,[23]). Thus, the cavity field 
will adiabatically follow the motion of its end-mirror, and 
so the bath creation and annihilation operators, b\ and 
6^, depend only on the index k (not explicitly on time). 

It is convenient to move to an interaction picture, in 
which the effect of the first two terms in equation ((T|) 
are factored into the energy basis states. In this rotated 
basis, the Hamiltonian takes the form 

Hi{t) =^2gk(t) &~b\ cxp J [uj k (r) ~ uj ]dT 



(2) 



Since the Hamiltonian commutes with N = a + a~ + b\.b k , 
the total number of energy quanta is a constant of the 
motion. We choose the atom to be in its excited state at 
t = and thus the state- vector at any later time is simply 
a superposition of those energy eigenstates corresponding 
to a single energy quantum: 



\Mt)) =c(f)|l><8>|...0...) 

+ 5^c*(*)|0>®|...lfc. 



(3) 



Using the Schrodinger equation, the equation of motion 
of the state-vector coefficients is found to be 



. dCgjt) 

1 dt 



.dc k {t) 
1 dt 



: YV (i)exp i-i [ 
u V Jo 



Wfe(r) - aj ]dT c k (t) 



(4) 



-g k (i)exp \ i J [w fc (r) - u a ]dT^ c a (t). (5) 



Since we are interested in controlling the atomic state, 
it is convenient to eliminate the environmental variables, 
Cfc(t), from equations ([4]) and ([5]) giving the following 
intcgro-diffcrcntial equation for the atomic state: 



f Q K{t,t')c a {t') dt', 



(6) 



with 



K{t,t') =J29k{t')g k (t)exp (-i J [u k (r) - <J ]dT^ . 

(7) 

We note that the atomic dynamics are determined solely 
by the behaviour of the kernel, K(t, t'), which in turn is 
dependent only on the properties of the atom's surround- 
ings: The principal idea underlying the use of dynamic 
environments is that by manually altering the reservoir 
properties, gk(t) and cjfc(i), it is possible to control the 
evolution of the atomic state. In other words, the in- 
evitable coupling between system and surroundings can 
be used to our advantage, since changes in the reservoir 
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have an effect on the atom, and this can be used to shape 
the atomic state. 

As a preliminary observation, we note that imposing 
a dynamic structure on the reservoir can only give rise 
to observable effects if the atom-field coupling is struc- 
tured, by which we mean that the product pkg 2 , varies 
significantly with frequency, where pk is the reservoir 
density of states. To see this, we note that in the con- 
verse situation of a completely structureless reservoir, 
with Pk9k{t) 2 = Po3o (constant), then whatever the mod- 
ulation function f(t) we choose for the mode frequencies 

El 



w fc (t) =wfc(0) + /(t), 
the kernel, (JT]), reduces to 



(8) 



Wfc(r) - LO ]d,T 



-9l 



2^exp 

k 



x exp 



(-*(w*(0)-wo)(t-t / )) 



f(r)dr 



i2-Kp gl5{t - 1'). 



(9) 



The population of the atomic excited state thus decays 
exponentially at exactly the rate predicted by Fermi's 
golden rule for a static reservoir: 



dCg(t) 

dt 



npo9o c a(t)- 



(10) 



Equation |[TU]1 holds, regardless of the frequency modula- 
tion f(t) that is imposed on the reservoir modes. There- 
fore, even though the atom is resonant with different 
field-modes at different times during its evolution, the 
overall effect of the reservoir is unchanged by modulating 
the mode frequencies alone. Dynamically modulating the 
microscopic properties of a reservoir can only have an ob- 
servable effect if the system-reservoir coupling is a struc- 
tured function of frequency. However, we note that the 
requirement of a structured reservoir is not synonymous 
with strong coupling between the system and its envi- 
ronment, since the distinction between strong coupling, 
reservoir structure and non-Markovian dynamics is more 
subtle for dynamic reservoirs than the usual static case. 
(For further details see [H[.) 



III. DYNAMIC RESERVOIR MODEL 

Having established that the use of dynamic environ- 
ments requires the system-reservoir coupling to be spec- 
trally dependent, we now turn to the specific form that 
this dependence will take. In order to isolate the novel ef- 
fects associated with dynamically structured reservoirs, 
we choose to investigate a model for which the macro- 
scopic properties of the reservoir remain fixed. (See figure 



Q]) This means that the individual coupling strengths, 
gk(t) must vary in time, to match the local reservoir 
structure at the changing frequency, u>k(t). In this way 
the envelope of the coupling strengths can be kept static. 
Put another way, the microscopic atom-mode couplings 
conspire in such a manner that the reservoir structure 
function Pk\9k{t)\ 2 must be expressible as a pure func- 
tion of u)k{t) only, with no explicit dependence on the 
time t. For a simple cavity model, we assume that the 



Envelope remains static 
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FIG. 1: This figure corresponds to the model studied in sec- 
tion [ITI] and shows reservoir mode couplings gk(t) (see equa- 
tion ([2])) at two different times ti and ti. The envelope of the 
couplings remains the same as the reservoir mode frequencies 
vary in time. Thus the coupling of the atom to an individ- 
ual mode changes in time. The thick line indicates the same 
mode k at two different times, ti and t2- Since at time t-z 
this mode has moved to a different frequency under the static 
envelope, the coupling has changed accordingly. 

single time reservoir structure function pk\gk{t)\ 2 is a 
Lorcntzian, with width 7, centred on the atomic tran- 
sition frequency luq, i.e. 



Pk\9k{t)\ 2 = 



D 2 "f/n 



7 2 + (w fc (f) - w ) 2 



(11) 



For a static bath, where Pk\9k{t)\ 2 is not a function of 
time, the Lorcntzian (1111) is a common choice for the 
reservoir structure [28l l29l| and the resulting dynamics 
have been well explored (see, for example, [3fJ, |31(). In 
the weak coupling limit, the Lorcntzian reservoir struc- 
ture ensures exponential decay of atomic population. 

Equation ([TTJ describes both the time-dependent cou- 
pling of a single bath mode, u>k(t), and also the instanta- 
neous coupling between the atom and the whole bath of 
modes at time t. Since this structure has no explicit de- 
pendence on t (i.e. only the implicit time-dependence 
contained in LOk(t)) the envelope of these couplings is 
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fixed in time; as in the static case, the atom always has 
some resonant modes with which it can exchange a pho- 
ton but these are now different bath modes at different 
times. We note that the weight, D, of the Lorentzian, 
which is defined through the relation 



(12) 



also remains static in time with this choice of couplings. 

Of course, the kernel, (O does not feature Pk\9k{t)\ 2 : 
but the two-time product Pk9k(t)9k{t')- However, as al- 
ready mentioned, gk(t) may be chosen to be real. There- 
fore the two-time product follows immediately as 



Pk9k{t)gk(t') 



ttVT^+W*) ~^o) 2 ) 
1 

X vV + Mi')-"^) 5 )' 



(13) 



To proceed with the analysis, a specific form for the time 
dependence of the modes must also be chosen. We as- 
sume that all the reservoir modes experience an identical 
frequency-modulation, which is periodic in time, with pe- 
riod SI 



w*(f) =wfe(0) + /(*). 



(14) 



As well as the modulation frequency, Q, it will also prove 
helpful to introduce the modulation depth 
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(fu 



/min) j 



(15) 



in order to characterise the amplitude of the frequency 
manipulation. We also note that the specific form f(t) = 
dsm(flt) will be considered in detail in section [TV] 

By choosing f(t) to be periodic, we ensure that the 
dynamically acquired phase-term arising from the reser- 
voir manipulation is also periodic, and thus has a discrete 
Fourier-decomposition: 



cxp 



f(r)dr 



E 



F n e 



-inVtt 



(16) 



with 



F„ 



cxp ( -i I f{r)dr I e~ inUt dt. (17) 



For this specific frequency manipulation wc note that the 
Hamiltonian j2|) may be re-expressed as follows: 



fl>(*)=5>*(*) 



a S[exp ( i (wfe(O) - u )t + i f(r)d 



+ a + b k cxp ( -i (w fc (0) -uo)t-i J /(r')dr' 



= E { E^^W^'^e'^^^^^^+KfffcW^^^ 1 ^^^ ""^^. 



(18) 



n— — oo ^ k 



Written in this way, we see that the chosen dynamic reser- 
voir structure (figure [l} is effectively the same as an in- 
finite collection of independent reservoirs, indexed by n, 
which have static mode-frequencies, and time-dependent 
couplings, F n gk(t). (These are the terms in curly braces 
in equation (fTS)) .) This picture is illustrated in figure 
[2j The n th sideband peak in the reservoir structure has 
central frequency u>o + nfl. From the definition of the 
coupling strengths, (fTTj) . we see that the n th sideband 
has the largest coupling when 



f(t) + nfl «0. 



(19) 



Since the condition, (fH?)) is different for each sideband, 
the atom decays into different sidebands at different 
times. The total decay-rate is therefore determined by 
the decay-rates into each individual reservoir, as shown 
in the following section. 

A. Atomic decay rates 



Using equations (|13j) and (|14)) 
explicit expression for the kernel 



we can now write an 
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K(t, =*1 £ expHMO)-^-t')] . exp r f{r)d \ {2Q) 

npk k y(7 2 + K(o)-^ + /(i)) 2 ) ( 7 2 + K(0)-c + /(t')) 2 ) V 7t ' 7 

It is convenient now to move to a continuum limit for the mode frequencies. For economy, we choose to write the 
initial frequencies simply as Wfc(O) = u>, which gives: 

r,, a ^ 2 7 /"°° exp [— i (u - uj q ) (t - 1')] , / r\r^\ 

K(t,t') = / PL v =rfu • exp -« / f(r)dT . (21) 



coupling 




Range of f(t) for lowT 
« ► 



Range of f(t) for high r 

FIG. 2: By re- writing the Hamiltonian (TJJ in the form (JTSj , 
the model studied in this paper may be interpreted as describ- 
ing the decay into a collection of reservoir sidebands, evenly 
spaced by angular frequency Q. In this picture, all reservoir 
mode frequencies are static in time, but the coupling profile 
of each bath is modulated according to (|13|) (vertical arrows) . 
Each bath passes through resonance at a different time, spec- 
ified by equation (|19p . As discussed in below equations (|3U 
I33[) . the atomic decay-rate, Too, can be increased or decreased 
by tuning the extremal values of f(t) towards or away from 
the sideband frequencies, \f(t)\ fa nQ. 

The kernel contains information about the 'memory' 
of the atomic dynamics, in the sense of equation (|6|). 
As a special case, we note that when the kernel decays 
sharply away from t' ~ t, memory effects are absent and 
the Markov approximation can be applied, to give: 



dCgjt) 

at 



'-Ca(t), 



with 



r- 



= lim <^ -5R 

t^oo 1 t 



t rtt 

dt\ I dt-2 
Jo 



K(t u t 2 ) 



(22) 



(23) 



For the Markov approximation to hold well, the Ker- 
nel must decay on a time-scale which is very short com- 
pared to the typical time-scale of atomic dynamics. For 
static reservoirs this requires that the system-reservoir 
coupling is weak (i.e. D <C 7). However, for dynamic 



reservoirs the atomic behaviour can be Markovian, even 
when the atom-reservoir coupling is strong, as long as the 
two-time product in the denominator of K{t\, t?) decays 
sufficiently fast. Thus, for dynamically manipulated en- 
vironments, the atomic behaviour can be Markovian even 
in the strong-coupling regime j4j|. The results that fol- 
low in this paper therefore hold when either of the two 
below conditions are satisfied: 

1. Weak-coupling: fl«7 

2. Fast rate of change of reservoir mode- 
frequencies: dO<C7 2 ,.D 2 . 

Our task is now to solve for the atomic decay-rate, 
Too 1 and to see how this can be controlled by altering the 
properties of the modulation- function f(t) which charac- 
terises the environmental manipulation. To this end, we 
note that (|23[) can be written in the compact form: 



D 2 7 



t— too 7rt 



lim 



J -oo \Jo 



x ( I dh P{ojM) 



(24) 



with 



P(u>,t) = 



exp 



-ifo /(">") dr 



^ 7 2 + (w-^ + /(t)) 2 
x exp [— i (uj — ujq) t] . 



(25) 



The mathematical steps involved in proceeding from 
equation (|2"3"|) to equation (TJH) include a change of lim- 
its, a re-assignation of dummy indices and use of the fact 
that Too is defined to be real. The term J* 1 P{u>,ti) dt\ 
in equation ([24]) is thus only significant if the two oscillat- 
ing factors which make up P(u>, t) interfere constructively 
over many cycles. Due to the periodicity of f(t), the 
term in curly brackets in equation (|25[) oscillates with pe- 
riod 27r/fi. The condition for the other oscillating term, 
exp [—i (u> — loq) t] to remain in phase with this term is 
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u>o fisnft with 



0,±1,±2,. 



(26) 
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Reservoir modes with frequencies which do not satisfy 
this condition do not contribute significantly to and 
therefore we should expect that they do not become oc- 
cupied during the atomic decay process. This is indeed 
found to be the CclSG, clS outlined in subsection Iff! Bl 
With the above interference argument in mind and an- 
ticipating that only those frequencies satisfying (|26[) will 
contribute significantly to T^. we choose to re- write 
equation (|24|) as follows: 




t 2 
P(loq + nil, t\) dt\ 



(i 



(27) 



The proportionality constant in (|27p is independent of 
the index n. In order to fix the normalisation of r M , we 
note that the function P(tuo+n£l, t) oscillates with period 
27r/0, and so we may replace the long-time average by 
the average over a single period: 



It is worth mentioning that although we chose to move to 
a continuum picture from the modes, we have now arrived 
back at a discrete sum over reservoir variables. However, 
now the index (n) runs over reservoir sidebands (which 
are collectively defined), rather than individual modes 
(k) with which we began. 

It is clear from equations (|30|) and ([31]) that in or- 
der to reduce the total decay-rate, the individual decay- 
channels must be inhibited, and that this can be achieved 
for the n th channel by making sure the following two con- 
ditions are satisfied: 

1. g( n '(t) is small while ui^ n '(t) is slowly- varying, 

2. oj( n \t) varies rapidly while g^ n \t) is large. 

This means that the atomic decay-rate can be inhibited 
if we choose the modulation function in such a way that 
\f(t)/fl\ is far from all integer-values at turning-points of 
f(t) (see figure[2]), and also f2 3> 7. 



Z? 2 fi 7 ^ 
I 00 oc - „ 2 „ 



2tt 2 



n— — 00 



P(lUq + 7lf2, t\) dtj 



Finally, we take the dual limit: 

f2 — > 00 and d — > 0, 



(28) 



(29) 



which corresponds to a static reservoir, and must there- 
fore return the corresponding (normalised) static decay- 
rate, = 2D 2 /j. This gives: 



r « V r 

1 00 ~ 7 ^ n-) 

n— — oo 

with the following three definitions: 



(30) 



Fn = 2tt 



g {n> {t)exp -i / w W(r) -uj 
\ Jo 



dr dt 



w w (t) = wB + «n + /(t) 



D 2 7/7T 



7 2 + ( w W(i) -cu y 



(31) 
(32) 

(33) 



B. Reservoir occupation-spectrum 

In order to better understand the mechanism by which 
the atomic decay rate takes takes the form in equation 
(p0| . it is useful to consider the final occupation spec- 
trum of the reservoir modes. We define this spectrum 
by considering those modes k$ (kg £ {k}) which have 
frequencies that lie within 8u> of an initial frequency, u>. 
That is, we let [H [H 



S(u>) 



lim 



i? Mf) ' 2 } 



(34) 



In the continuum limit the number of modes kg in the 
sum is approximately 6uip(u>), and since the cu are ex- 
pected to vary smoothly with k in this limit, we finally 
let 

> lim \p(w)\c ks (t)\ 2 }, (35) 

which applies to a representative kg. Equation (|35[) will 
serve as an operational definition of the bath spectrum. 
From equation ([5]) we find the corresponding solution for 
Cfc(t), in integral form: 



lim \c k {t)\ 



1 I g k (t)ex.p 



cjfe(r) - uj Q + 



dr 



dt, 



/o V 2 

which, together with the definitions (fT4|) and (fT3|) and the Fourier decomposition (fT7|) gives 

S(u) - 



POO 

/ dt ^2 F, 

«/ — _ „ 



exp [i (uj — luq - mft + ^f 2 -) t] 



m— — oo 



^^ + {uj-uj + }{t)f 



(36) 



(37) 
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The interference arguments used in deriving an cxprcs- occupation spectrum can be written as a discrete collec- 
sion for can be applied again here, and so the final 
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FIG. 3: Final reservoir occupation spectrum for the decay 
process studied in section IIII Al The vertical axis shows the 
discretised form of S(lu), i.e. limt^oo {p/c|cfc(t)| 2 }, calculated 
via a numerical solution to equations @ and ((SJ). The ex- 
ample shown uses D = 7, SI = 2O7, d — 687. (a) The final 
occupation spectrum is a discrete collection of peaks, centred 
at wo + nQ,. This is due to the dynamical resonance effects 
discussed in section ITTT1 The area bounded by a peak indexed 
by n is r^/Too. (b) The shape of each peak is approximately 
Lorentzian with a half- width Foo/2, as predicted in section 

nxmi 



tion of peaks at frequencies: 



LO — UJq 



inSl 



with 



0,±1,±2,. 



(38) 



Figure[3l shows the final occupation spectrum, plotted 
using results from a numerical simulation of equations ^ 
and (O for the specific case of a sinusoidal modulation 
(studied further in section HV]) . which shows that the fi- 
nal spectrum is indeed a discrete collection of Lorentzian 
peaks, as predicted. The effect of dynamically manipu- 
lating the individual frequencies and coupling constants 
is to create a sideband structure for the reservoir. Each 



sideband, indexed by the integer n, is centred at the fre- 
quency given by equation (f3"8"|) . 

We write the weight of the n th peak as 

ru> +(n+l/2)Q 

S n = S(u) du, (39) 

L +(ii-i/2)n 

and (largely following the arguments in section IIII Ap it 
is reasonably straightforward to show that 



O 2 



27rr c 

no 



3 W (t)exp ( -i J uj {n > (t) - luo 



dr dt 



(40) 
(41) 



IV. EXAMPLE: SINUSOIDAL MODULATION 
OF RESERVOIR MODE FREQUENCIES 

A. Decay with sinusoidal oscillation 

The analysis of the previous section is fairly general, in 
that the sole requirement placed on the reservoir manip- 
ulation is that the modulation function f(t) is periodic 
in time. In order to present a concrete test of the pre- 
dicted decay-rates and reservoir occupation spectrum in 
the current section, we apply our results to the simplest 
periodic case - namely a sinusoidal frequency manipula- 
tion, with modulation depth d and modulation frequency 
O 



/(f) =dsin(ftt). 



(42) 



For this case, the Fourier decomposition of the phase- 
term is : 



exp [id (cos(nt) - 1) /SI] = e - ld ' n *" J n (d/Sl)e inQt , 

n— — 00 

(43) 

where J„ is the n th -oidei Bcsscl function of the first kind. 



The decay-rate into the n th sideband given in equation 
(pTTj) thus takes the form: 



D 2 Sl 2 j 
2tt 2 



J2 * m J,n(d/Sl) 



m— — 00 



-i(n — m)Qt 



yV + (nSl + dsm{Slt))' 



dt 



(44) 



r 



The total atomic decay-rate, = J2 n ^« ^ s plotted 
in figure [5] and compared with a direct numerical solu- 



tion of equations (0| and |5]) for a range of modulation 
depths and fixed modulation frequency. As explained in 
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FIG. 4: The result of equation (|30[) is plotted as a red line 
and compared to numerically extracted decay rates, shown 
as blue crosses. The numerically extracted decay-rates were 
calculated by solving equations @ and © and fitting an 
exponentially decaying curve to the resulting |c a (£)| 2 . The 
modulation frequency Q — 2O7. 




all other frequency-scales in the system (except for the 
free-evolution) : 



iv > tt > d,7,£>. 



(45) 



In this limit, the denominator of every term is small for 
n 0. Also, the sum over m is dominated by the m = 
term, due to the factor J m (d/Q), and so we only need 
to consider the n = m = term. This can be tackled 
analytically: 



2tt 2 



2D 2 



Md/Q) 



1 



dt 



2 7 j (d/ft)x; * 



Y + <P sin 2 (m) 
2 ~\ 



7T (7 2 + <P) 



(46) 



where fC represents the complete elliptic integral of the 
first kind |33j |. In this limit, is a monotonically de- 
creasing function of d/7 and the inhibition of atomic de- 
cay can be quite dramatic, as illustrated in the following 
section. 



V. PHYSICAL REALISATION AND 
OBSERVABLE EFFECTS 

A. Double-cavity system 



FIG. 5: Final occupation peak-weights, S n as a function of 
the modulation-depth, d for the first four peaks. The figures 
(a)-(d) correspond to n = 0, 1, 2, 3 respectively. The n th peak 
weight is predicted to be r n /Too, with Y n given by (|31|l . As 
discussed in subsection lHI Al S„ ~ for d < nil and S n peaks 
at d ~ nf2. 

the text of section IIII A[ is a peaked function of the 
modulation depth, with a maximum every time d is an 
integer- multiple of £1. 

The final reservoir occupation spectrum, S(w), calcu- 
lated numerically from a direct simulation of equations 
(HI) and is plotted in figure EH As predicted, this spec- 
trum is a collection of Lorentzian peaks centred at angu- 
lar frequencies u m ljq + nfl. 

In order to test the predictions for the weight, S n = 
Tn/Foo, of each peak, figure O shows the analytic result 
of equation (f4"Tj) together with the numerically calculated 
peak weights defined in equation (|39p , for the same model 
parameters as figure [U 



The reservoir manipulation studied in this paper is 
quite specific in the sense that the microscopic reservoir 
modes have time-dependent frequencies and couplings to 
the atom, but together, they must conspire to keep the 
macroscopic reservoir structure fixed in time. This pre- 
cise form was deliberately chosen in order to isolate the 
new effects arising here from other prev iously studied 
types of engineered reservoir, such as |l8l . l3~il . 13a . [3(| 

In a recent paper, we suggested one way of physically 
implementing such a dynamic environment [23| ; the atom 
is placed inside a static cavity, which itself is enclosed 
in a larger cavity with moving boundaries, as in figure 
(S](a). In this way, the envelope of the reservoir structure 
is determined by the properties of the inner cavity, and 
so remains fixed. On the other hand, the frequencies of 
individual modes are fixed by satisfying a node condition 
on the outer mirrors, and so these can be modulated by 
moving the right-hand mirror. 

B. Comparison with a variable-detuning model 



B. Decay in the limit of ultra- fast reservoir 
modulation 

The simplest analytic expression for Too is obtained in 
the limit that the modulation frequency dominates over 



It was noted in section [TV] (e.g. figured]) that dynami- 
cally manipulating the reservoir structure in the way de- 
scribed here can result in an inhibition of the atomic 
decay rate, despite the fact that the atom and reservoir 
are weakly coupled, and always remain on resonance dur- 
ing the decay. In this subsection, we choose to highlight 
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Moving Mirror 
« - > 




> u ' Moving Mirror (0 

< ► 




L,(t) 



OflT-re sonant 
laser field 



FIG. 6: (a) To manipulate the microscopic mode-structure 
of the reservoir (double-cavity) , the outer-cavity mirrors are 
moved, which affects the mode frequencies as described by 
(|14l) . Since the inner-cavity mirrors are fixed, the macroscopic 
properties of the reservoir remain invariant which recovers the 
Hamiltonian of equation {TJ. (b) Modulating the length of 
the cavity itself acts to dynamically detune the cavity reso- 
nance from the atomic transition, resulting in the interaction 
Hamiltonian (|4T[) . (c) Alternatively, variable detuning be- 
tween the atom and reservoir may be achieved by inducing a 
time-dependent Stark-shift in the atomic transition frequency 
[37| . This method also results in the Hamiltonian given in 

SB- 



that the extent of this inhibition is far greater than might 
first be expected. To this end, we choose to contrast the 
situation studied in sections HIllIVI with with the case of 
a variable atom-reservoir detuning, i.e. 

H'i{t) =y^gfc <J~6j.cxp^i J [uJk(r) - uj ]dT 



+ & + b k exp ( -i J [w fc (T') - cj ] dr' 



(47) 



Modulating either one of the atomic transition frequency 
(figure Efc)) or the central frequency of the reservoir (fig- 
ure EJb)) would give rise to the interaction Hamiltonian 
in equation (|47p . In either case, the coupling strengths gu 
are static and do not compensate for the changing mode- 
frequencies and this has the effect of dynamically mod- 
ulating the macroscopic atom-reservoir detuning. Such 
variable-detunin g m odels have been well studied (see for 
example (2ll 123, 123, l38j]) and in the weak-coupling limit 
the atomic decay-rate is given by j22| : 



r' 



r' 



2D 2 



-j n (d/n) 2 



7 



7 2 + n 2 Q 2 



(48) 
(49) 



In the ultrafast modulation limit considered in section 
HVB[ equation (]49fl simplifies to 



r' 



2D 2 

1 



(50) 



Somewhat surprisingly, in the limit of ultrafast frequency 
modulation, the decay-rate T'^ is always greater than the 
equivalent decay rate given in equation (|46|) for the 
same frequency modulation and depth. Therefore a more 
effective suppression of the atomic decay-rate can be 
achieved by keeping the atom and reservoir on-resonancc 
rather than dynamically changing the detuning. 

The approximate analytic results given in equation 
PS]) and equation (|5D)l are compared in figure [7] together 
with the full decay-rates given in equations (|46|) and (|49|) . 




FIG. 7: Analytic results, (J44J) and {49} (crosses), together 
with the leading-order terms, (|46|l and (|50[) (solid lines). The 
modulation frequency SI — 2O7. In this limit, the reservoir 
manipulation studied in our paper (lower curve) is more ef- 
ficient at suppressing dissipative effects than an oscillatory 
detuning (upper curve). 

We note also that the ratio between the decay rate 
derived in this paper (section IIV Bjl and the decay rate 
arising in a variable-detuning model may be written as 



r' 



K 



(its 7 ) 



TT 2 (l + X 2 ) 



(51) 



with x = d/7. This is a monotonically decreasing func- 
tion of x, and always less than 1. Providing that the 
modulation frequency is high enough, it is possible to 
achieve high values of x whilst still satisfying ft 3> d. 



C. Size of effects 

Experimental verification of the predictions made 
above may be achieved with currently available tech- 
nology and without great difficulties. First we take 
an optical atomic transition with angular frequency 
lu /2it ~ 3.5 x 10 14 Hz and a 40 /im inner cavity for 
which 7/2?r - 4.1 MHz [H, SO, EH, whilst for the outer 
cavity we let L = 1cm. Assuming that the end-mirror is 
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driven sinusoidally by a piezoelectric actuator, with fre- 
quency 0/27T ~ 40MHz and an amplitude d ~ 0.4nm, 
the ratio f2 > d > 7 is achieved and the resulting decay 
rate, is four times lower than for a variable-detuning 
model. 

For microwave systems, we consider an atomic transi- 
tion with frequency lvo/2tt ~ 21GHz and a 2cm inner cav- 
ity for which r y/2n ~ 10 Hz Q, whilst for the outer cavity 
we let L ~ 20cm. In this case, a frequency il/2n ~ 1MHz 
and an amplitude d ~ lOnm, gives Sl/d = d/7 = 1000 so 
the heirachy given in equation (|45[) can be made to hold 
very strongly. The resulting inhibition of the decay-rate 
(equation ([BTj)) is ^ ~ 2.8 x 10 _5 rj x) ; decay into the 
cavity field mode is effectively 'switched-off' in this limit 

VI. CONCLUSIONS 

We have studied how the process of spontaneous emis- 
sion may be controlled by dynamically manipulating the 
microscopic structure of an atom's environment. While 
is it well known that changing the structure of an atom's 
surroundings (by introducing nearby conducting surfaces 
or a cavity, for example) can alter the process of spon- 
taneous emission, in this article we have considered a 
model for which the macroscopic reservoir-structure does 
not change with time. Instead, the individual electro- 
magnetic field-modes which make up the reservoir are 
assigned time-dependent frequencies and time-dependent 
couplings to the atom. In this way, the atom and reser- 
voir always remain on resonance but the atom interacts 
with different bath-modes at different times; by manip- 
ulating the environment in this manner it is possible to 
control the memory-kernel for the atomic state and hence 
to control the atomic decay-rate. 

In order for such a scheme to work, we must be able 
to alter the reservoir-structure on a timescale which is 
shorter than the static-bath correlation time. There- 
fore, it is essential that the system-bath couplings are 
frequency-dependent, as shown in section [TTJ Although 
essential, this frequency-dependence may be very weak, 
and in section HTTl we show that it is possible to control the 
atomic decay-rate even in the case of weak atom-reservoir 
coupling (as well as in the strong-coupling limit - see 
[23l]). In subsection IIII Al the atomic decay-rate is calcu- 



lated for the case of a periodic manipulation of all reser- 
voir mode- frequencies. A Floquet analysis is applied and 
after analysing the final occupation-spectrum (in subsec- 
tion llll B|) , we find that the atomic decay only occurs into 
a discrete collection of reservoir-sidebands, each of which 
satisfies a dynamic resonance condition with the atom. 
The total atomic decay-rate is a structured function of 
the modulation-depth and modulation-frequency and ex- 
hibits a peak whenever one of the reservoir sidebands is 
resonant with the atom for a significant fraction of each 
modulation cycle. 

Section ITVl treats the specific case of a sinusoidal mod- 
ulation of all reservoir mode-frequencies and the general 
expressions derived in section imi are found to be in excel- 
lent agreement with numerical simulations. Of particular 
interest is the limit of ultrafast modulation, studied in 
subsection IIV Bl for which compact analytic expressions 
are derived. 

In section [V] the observable effects of the model are 
considered. In order to emphasise how strong the sup- 
pression of decay may be in the limit of ultrafast mod- 
ulation, a second model is introduced in subsection IV Bl 
for which the atom-reservoir detuning is modulated sinu- 
soidally. Upon first inspection, it may appear that this 
second reservoir manipulation should give rise to a more 
dramatic inhibition of the atomic decay, since in this case 
the atom and reservoir spend a large proportion of each 
cycle away from resonance. However, a detailed analysis 
shows that the dynamical suppression of decay studied 
in sections Mil and IIVI is a far more potent method for 
inhibiting dissipation than a straightforward atom-bath 
detuning. Physical realisations for both of these types 
of dynamic environment are proposed in subsection IV Al 
(also, more details are given in [23|| ) . Finally, in subsec- 
tion |VCl we show that it should be possible to observe 
all of the effects predicted above using currently available 
technology. 
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